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Abstract 

Mean flow effects are discussed for two different pattern-forming systems: Rayleigh- 
Benard convection and Faraday instability in viscous fluid. In both systems spirals 
are observed in certain parameter regions. In the Rayleigh-Benard convection, the 
spiral core instability and subsequent generation of up- and downflow hexagons are 
shown to occur due to the mean flow generated by the curved rolls near the core. 
In the Faraday instability, the mean flow which is generated by rapidly decaying 
surface waves near the wall, causes wavenumber frustration which leads to a rigid- 
body spiral rotation. In both cases we use phenomenological Swift-Hohenberg-type 
equations for the order parameter coupled to a large-scale mean flow. Numerical 
simulations are compared to recently reported experimental results. 



1 Introduction 



Nonlinear evolution of cellular patterns in extended systems near the threshold 
of a primary instability are usually studied using amplitude equations (e.g., 
Newell- Whitehead-Segel equation for roll patterns in isotropic systemsfl]). 
These amplitude equations can be derived from the first principles by pertur- 
bative expansion using the supercriticality e as a small parameter. However, 
the expansion is based on a pre-chosen orientation of the cellular pattern and 
thus do not preserve the full spatial symmetry of the physical system. There- 
fore this description fails for more complicated patterns such as coexisting do- 
mains with varying orientation of rolls, targets, spirals, etc. To overcome this 
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difficulty, phenomenological models similar to the so-called Swift-Hohenberg 
equation (SHE), 



are often employed (see [2] and references therein). This model correctly de- 
scribes the linear properties of the system close to the threshold and yet 
preserves the rotational symmetry of the system. The disadvantage of this 
approach is that the nonlinear term in this equation is usually added ad hoc 
in order to provide a saturation of the primary instability and cannot be sys- 
tematically derived from underlying first principles. Thus, while some of the 
features (including some secondary instabilities, domain coarsening, etc.) of 
the pattern formation can be modeled based on this equation, some other 
essential features of the dynamics may be missing. In many examples these 
features are related to the coupling between the cellular mode and so-called 
zero mode slowly varying both in space and time. It can be a concentration 
field in binary-fluid convection [3], population inversion in extended lasers[4], 
density field in granular media [5], mean flow in Rayleigh-Benard convection [6- 
8]. This problem can be rectified by expanding the basic model (1) so that it 
includes the interaction with the zero mode. It allows to reproduce many subtle 
features of the dynamics of cellular patterns within a simple phenomenological 
description. 

In this paper we illustrate these general ideas by two examples of spiral dy- 
namics in Rayleigh-Benard convection and oscillating fluid layer. In both cases 
the spiral rotation is intimately related to the interaction with mean flow, self- 
generated by the cellular patterns in case of convection, and wall-generated in 
case of Faraday system. We show that recently observed spiral core instabil- 
ity[9,10] and spiral-hexagon transition[10] in Rayleigh-Benard convection in 
Boussinesq fluid can be understood within the model for the order parameter 
coupled with the mean flow [11]. We also demonstrate that spiral rotation in 
Faraday experiment can be explained by the wavenumber frustration which 
is caused by the near-wall radial mean flow generated by decaying capillary 
waves [12]. 



2 Spirals in Rayleigh-Benard convection 

Recently discovered spiral-defect chaos (SDC) in large aspect ratio Rayleigh- 
Benard convection system at small Prandtl number[13] was successfully repro- 
duced both within Navier-Stokes equations[14] and also using SHE coupled to 
the equation for the mean flow generated by the curved rolls[15], 
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Here ip is the order parameter, u the horizontal velocity field of the large-scale 
flow, and ft the vertical component of the vorticity. The control parameter 
e represents the reduced Rayleigh number, while a characterizes the Prandtl 
number of the fluid. The parameter g allows to more accurately reproduce the 
stability properties of convection patterns [8], and g m characterizes the cou- 
pling strength between the order parameter ip and the vorticity f2. The phe- 
nomenological parameter c is introduced to describe the local dissipation of the 
vorticity (e.g. due to friction at the bottom of the convection cell) [8,2]. Thus, 
Eq. (2) describes the dynamics of the order parameter if), while Eq. (3), using 
the definition of the vorticity ( 4), represents the coupling of the large-scale 
flow field u and the order parameter. For g — 1 and g m = Eqs. (2)- (4) reduce 
to the Swift-Hohenberg equation (1). We solved Eqs. (2)-(4) numerically in a 
domain of 256 x 256 collocation points using a pseudo-spectral method based 
on the Fast Fourier Transform. The physical domain size was typically re- 
stricted to 150 x 150. Circular boundary conditions were enforced by ramping 
e towards negative values at distances from the center r > -Rmax = 55. 

Essential features of SDC include spontaneous spiral creation, quasi-stationary 
spiral rotation, spiral core instability, and eventual spiral destruction by other 
spirals. Cross and Tu[16] showed that persistent spiral rotation is caused by 
the wavenumber frustration between the values selected by the spiral core and 
by the environment (most typically, roll dislocations). It is noteworthy that 
large vorticity generated in the spiral core plays only the secondary role in 
this mechanism since the asymptotic wavenumber selected by the core does 
not depend on the Prandtl number and coupling to the mean flow (parameter 
g m in Eq.(3)). However, vorticity plays a major role in wavenumber selection 
by the roll dislocation. It is easy to see that the dislocation in the roll pattern 
generates a vortex pair which drives the dislocation towards the half-plane 
with greater wavenumber (see Fig.l where the order parameter for individual 
spiral terminated by a dislocation is shown together with corresponding mean 
flow vorticity). At small g m e the dislocation climbing speed due to vorticity is 
v oc g m e. The dislocation is stationary when this mean-flow drift is balanced 
by the ordinary climbing velocity v c oc (1 — gd) 3 / 2 [17], therefore the selected 
wavenumber q d0 oc 1 — a(g m e) 2 ^ 3 (a is a constant depending on other parame- 
ters g, a, c). In a stationary rotating spiral, wavenumber q d near the dislocation 
is larger than q d0 so the dislocation moves around the core of the spiral with an 
angular velocity coinciding with the frequency of the core rotation. It is easy 
to see from the phase diffusion equation (see [16]) that the frequency of the 
spiral rotation is inversely proportional to the distance from the core to the 
dislocation. Detailed measurements of individual spiral rotation rate confirms 
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this prediction [18]. 

At large supercriticality, a novel instability of the spiral core was recently 
observed[9,10]. The core exhibits high-frequency oscillations (with a period 
of several vertical diffusion times as compared with a period of overall spiral 
rotation of a few hundred vertical diffusion times). Using phase approximation 
(see [2,16]) it is easy to see that this instability is caused by the large vorticity 
generated near the core of the spiral. Indeed, in the phase approximation, 
outside the core of the spiral the phase is described by 

e t + U- Ve = r(q)- 1 V- [qB(q)} (5) 



where r(q), B(q) are known functions of the wavenumber q = V9 (-B(l) = 
0, B(q d ) < for SHE), see [2]. For stationary rotating n-armed spiral, 6 = 
J q(r)dr + ncp — uit, and therefore local wavenumber q(r) obeys the following 
ODE, 

^ = -^-r(q)- 1 ^[qB(q)} } (6) 



where T is a total circulation of the mean flow around the spiral core. It is easy 
to see that for 1 — q^ <C 1 the wavenumber q(r) ~ 1 — Air — A 2 nT/r 2 where A 12 
are constants depending on r and B. The wavenumber slowly rises towards the 
core and then rapidly decreases in its near vicinity. Of course, this expression 
is only valid outside the core of the dislocation (r > 1), however it illustrates 
the tendency for the core vortex to locally unwind the spiral. If T gets large, 
the wavenumber near the core goes outside the stability balloon and the spiral 
core becomes unstable. Visually it appears as that the core rotates against the 
overall spiral rotation (see a series of snapshots in Fig. 2 and computer movie 
at |http: / /inls.ucsd.edu/ ~lev/conv/| ). This behavior is qualitatively similar to 



observed in recent experiment by Plapp and Bodenschatz[9]. The threshold 
of this core instability depends on both e and g m as well as the topological 
charge n, see Fig. 3. 

For yet larger values of supercriticality, a transition to hexagons in the core of 
the spiral occurs, see Fig. 4. Both up- and down-flow hexagons are generated 
simultaneously near the core of the spiral. The existence of stable hexagons in 
systems without quadratic nonlinearity such as Swift-Hohenberg which pre- 
serves the symmetry ip — > —i/j, was attributed to the spontaneous excitation of 
the zero mode by Dewell et al [19]. It should be noted that within the frame- 
work of pure variational SHE (1) the hexagonal state can only be metastable, 
so a domain of roll will always expand towards a domain of hexagons. In case 
of SHE coupled with mean flow, the zero mode is generated and stabilized near 
the core by the mean flow which unwinds the spiral and drives the wavenum- 
ber towards zero. This zero mode gives rise to a proliferation of hexagons near 
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the core of the spiral. This transition was observed experimentally in [10,11] 
(see Fig. 5). 



3 Spirals in a Faraday system 

Parametric instability of a flat fluid surface subject to vertical oscillations re- 
mains a popular experimental tool for various pattern formation phenomena. 
Squares, rolls, hexagons and higher-order quasi-patterns have been observed in 
different regions of parameters (viscosity, driving frequency, and depth of the 
fluid layer). Phase diagram of different patterns has been recently computed 
by Chen and Vihals based on the Navier-Stokes equations [20]. In the region 
of large viscosity and large h/X, straight rolls are a stable pattern. However, 
as in case of Rayleigh-Benard convection, stable rotating spirals have been 
recently observed in the same parameter region [2 1,12]. In this section we dis- 
cuss the origin of the spiral rotation within the order-parameter model for the 
parametric instability, 

^ = - v$ - (1 + ia) \ip\ 2 i/j + m(V 2 + - (u • V)^ (7) 

Here is a complex amplitude of surface oscillations at the parametric fre- 
quency uj (which is a half of the driving frequency), 7 is a forcing magnitude, 
k is the dispersion parameter, and u is the velocity of the mean flow. Linear 
terms in this equation can be derived from the dispersion relation for capil- 
lary waves under parametric excitation, expanded near u = u , q = 1. The 
nonlinear term cannot be derived rigorously, as in Eq.(2) for Rayleigh-Benard 
convection, and has been added ad hoc to account for the stabilization of the 
parametric instability. Imaginary part of the nonlinear coefficient a describes 
nonlinear frequency shift (see also [22]). 

The last term in the r.h.s. of Eq.(7) describes order parameter advection by 
the mean flow. Such mean flow was observed in experiments [12] near the walls 
of the cavity. This flow is caused by the momentum transfer from dissipating 
capillary waves at the driving frequency which are generated by oscillating 
side walls. This flow is directed off the walls near the surface and due to 
incompressibility returns back to the walls near the bottom Velocity u in (7) 
should be understood as an the average velocity over layer thickness weighted 
with the vertical structure of waves. Since surface waves decay towards the 
bottom, near-surface flow affects them stronger than near-bottom return flow 
and the net u is oriented towards the center of the cavity. 

Equation (7) with periodic boundary conditions was studied numerically using 
pseudo-spectral split-step method with 256 x 256 collocation points, domain 
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size d = 200 and integration time step 0.05. To simulate waves in circular 
cavity, we ramped linear dissipation outside the circle of radius r = 86, i.e. 
v = u Q ,r < r and v = z/ (l + k (r — r )), r > r , where k varied between 0.5 
and 1.0. We assumed that the flow had radial direction and was azimuthally 
symmetric, u = u(r)r. We used the following profile for flow velocity u(r) : 

u(r) = w exp[£(r - r )] (8) 

For 7 > v trivial state — is unstable with respect to perturbations 
with wavenumbers near 1. Numerical simulations show that at the nonlin- 
ear stage, these perturbations give rise to various cellular patterns, including 
plane waves, targets and spirals. Without mean flow term (u = 0), these 
patterns remain stationary even when nonlinear coefficient in (7) is complex. 
Nonlinear frequency shift oc a only leads to deviation of the selected wavenum- 
ber from q = 1. [In systems with ordinary (non-parametric) pattern- forming 
instabilities non-potential effects usually lead to wave propagation.] However, 
when near- wall flow (8) is introduced in (7), standing waves comprising tar- 
gets and spirals begin to drift slowly toward the center. Figure 6 shows rolls 
velocity as a function of uq and £. Note that £ ^> Tq 1 so the flow is absent in 
the bulk, still rolls are moving throughout the integration domain. The phase 
velocity of rolls grows linearly with uq as it should be expected. 

This phenomenology can be understood in terms of the phase diffusion equa- 
tion similar to Eq. (5) for Rayleigh-Benard system. The important difference is 
that here we have standing waves, and therefore two phases 9 + and should 
be introduced. However, for large parametric forcing, the sum of phases 
is slaved to the phase difference = 6+ — 0-. The equation for the latter reads 

d& + U • V0 = r-\q)V ■ [qB(q)}, (9) 

where q = |VG and r{q),B{q) are functions of the wavenumber q. For the 
spiral, 9± = ±(/ qdr + m0) — u±t. Multiplying (9) by r(q) and integrating 
from r = to -Rmax we obtain 

lu s J r(q)dr = - J U(r)qr(q)dr + qB(q)\ r=Rmax , (10) 

where uj s = — u + is a frequency of spiral rotation. For q w 1 the last 
term in (10) is small and the frequency of spiral rotation can be estimated 
as uj s oc JU(r)dr Mo£ _1 - This relation is in an agreement with numerical 
simulations (see Figure 6). 

The mechanism of spiral rotation which we described is based on the assump- 
tion that the radial flow is produced by the rapidly decaying waves near the 
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wall generated by the oscillating meniscus. If this hypothesis is correct, the 
phase velocity of waves should increase linearly with the magnitude of driving. 
It also should depend sensitively on the wall profile. Indeed, experiment [12] 
confirmed this prediction (see Figure 7). 



4 Conclusions 



In this paper we briefly considered two different pattern-forming systems 
within the framework of order-parameter models coupled with mean flow. 
It turned out that both for Raylerigh-Benard convection and Faraday system 
the mean flow plays a crucial role for spiral dynamics. In the former, vortex 
pairs generated by roll dislocations, move the selected wavenumber away from 
the value selected by the spiral core and lead to wavenumber frustration. The 
frustration in turn leads to the persistent spiral rotation. The strong vortex 
which is generated in the core of the spiral plays a minor role in the overall 
spiral rotation, however it leads to the local spiral unwinding and eventually 
to the spiral core instability. If the wavenumber near the core reaches zero due 
to local spiral unwinding, zero mode is generated and leads to proliferation of 
up- and down-flow hexagons. 

In the Faraday system, stable rotating spirals have been recently observed. 
We showed that rotation of these spirals is caused by the near-wall radial flow 
which leads to the wavenumber frustration in the bulk. This near-wall flow is 
generated by rapidly decaying meniscus waves at the driving frequency. 

The results presented in this paper have been obtained in collaboration with 
I.Aranson, M.Assenheimer, and V.Steinberg (Rayleigh-Benard convection) and 
S.V.Kiyashko, L.N.Korzinov, and M.I.Rabinovich (Faraday system). 
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Fig. 1. Snapshot of the order parameter tp (left) and vorticity O (right) for 
a stationary rotating spiral terminated by a dislocation for Eqs.(2)- (4) with 
e = 0.7, g = 1, g m = 50, c 2 = 2, a = 1. 
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Fig. 3. Stability diagram for one- (a) and two-armed (b) spirals for c 2 = 2, a = 1 and 
g = 1. As e increases, the spiral core becomes unstable at the dashed line, regaining 
stability at the solid line as e decreases. 
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Fig. 4. Hexagon nucleation at a spiral core obtained for Eqs.(2) -(4) with 
e = 1.9, (7 = 0.75, g m = 10, c 2 = 2 and a = 1. Snapshots taken at 
i = 10, 110, 470, 650, 340, 2350. 




Fig. 5. Hexagon nucleation at a spiral core in experiment with e = 3.19, a = 4.5 
and time delay between frames At = 3.6, 3.6, 22.7, 18.0, 10.7 t v . 
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Fig. 6. Phase velocity of waves as a function of near- wall flow magnitude no (a) and 
inverse scale £ (b) from numerical simulations of (7). Parameters of simulations: 
7 = l.o, v = 0.5, a = 0.0, k = 1.0, fc = 1, £ = 0.1. In a, £ = 0.1, in 6, u = 1.0. 
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Fig. 7. Velocity of wave drift (measured by nodes displacement) as a function of mag- 
nitude of vertical acceleration for two different profiles of side walls. Inset: Sketches 
of the vertical profiles of the side wall and the structure of the shear flow. 
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